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ABSTRACT 

Analyses are presented for compressible fluid flow across shaft face 
seals with face deformation. The solutions are obtained from ah approxi- 
mate integral analysis. The models, used in this analysis, can predict gas 
film seal behavior operating at subsonic or choked flow conditions. The 
flow regime can either be laminar or turbulent. Entrance losses can also 
be accounted for. When fluid inertia effects are negligible (quasi-fully- 
developed flow) and the sealing faces are slightly deformed, the following 
results are found for both laminar and turbulent flows: 1) The pressure 

profiles are independent of fluid properties; 2) The parallel film leakage 
equation can be used, provided a characteristic film thickness is used. 
However, fluid inertia effects were found to be very important for near- 
choked and choked flow conditions. For these conditions the pressure pro- 
files are dependent upon the fluid properties. Pressure profiles are pre- 
sented for both divergent and convergent seal faces under choked flow 
conditions. 

NOMENCLATURE 
2 2 

A cross-sectional area, in ; m 

C T velocity entrance loss coefficient 

D hydraulic diameter, 2h 

f mean Fanning friction factor, 


/ Pu 
/ 2 


h film thickness (gap), in; m 

2 2 

h , characteristic film thickness, (h,h~/h ) , in: m 

char 12m. 

i specific enthalpy 

M Mach number 

M mass flow, lbm/min; kg/sec 

2 

P pressure, psi; N/m 

sealing dam inner radius, in; m ' 

sealing dam outer radius, in; m 
|R gas constant, ft-lbf / (lbm) (R) ; J/(kg)(K) 

Re leakage flow Reynolds number in radial direction, pUh/y 

r radial direction coordinate 

r f 1 linear tilt factor, h^^h^ + otx)/2h m (h^ + ax)^ 

T temperature, F; K 

u mean velocity in r-direction or x-direction, ft/sec; m/sec 

W flow width, in; m 

x coordinate in pressure gradient direction 

z coordinate across film thickness 

a relative inclination angle of surfaces, m rad 

0 radial flow factor, W/r 

y specific-heat ratio 

2 2 

p absolute or dynamic viscosity, (lbf ) (sec) /in ; (N)(sec)/m 

2 4 3 

p density, (lbf ) (sec ) /in ; kg/m 

subscripts: 

char characteristic 

h based on film thickness 

w wetted surface 


2 


0 


sealed (reservoir) conditions 


1 entrance conditions 

2 exit conditions 

3 ambient sump conditions 
superscripts: 

* referenced to Mach one condition (critical flow) 

INTRODUCTION 

Shaft seals in advanced aircraft rotating machinery will operate at 
speeds, temperatures, and pressures higher than shaft seals currently used. 

An example is shaft seals for advanced aircraft turbine engines. Conven- 
tional face contact seals presently used in gas turbine engines are limited 

to sliding velocities of about 350 feet per second (110 m/sec), pressure 

2 

differentials of about 125 pounds per square inch (86 N/cm ), and gas tem- 
peratures of 800 F (700 K)(l); pressure and speed capabilities of circum- 
ferential seals are near that of face seals. Advanced engines, however, 

will require seals to operate to speeds of 500 feet per second (150 m/sec) 

2 

(2), pressures to 500 pounds per square inch (340 N/cm ) and temperatures 
to 1300 F (980 K) (3) . Because of these severe operating conditions, seal 
face deformation is very likely to occur. 

These deformations may be due to various distortions (thermal, cen- 
trifugal, pressure, etc.). Seal face distortions become more pronounced 
under severe operating conditions and are usually detrimental to seal per- 
formance. Hence, prediction of these face deformation effects on gas film 
seal performance is of paramount importance. 

For face seals operating under severe conditions, a positive face separa 
tion (no rubbing contact) will be required in order to achieve long life and 
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reliability. A successful method of maintaining positive seal face separation 
is to add self-acting lift pads, such as shrouded ' Rayleigh step bearings, to the 
conventional pressure balanced face seal (4 and 5) . This is illustrated 
in Fig. 1. The self-acting lift pads (gas bearings) have a desirable char- 
acteristic, a decrease in film thickness results in an increase in the 
opening or separating force. Thus, the pads give axial film stiffness to 
the seal so that the stationary nosepiece will dynamically track with the 
rotating seal seat. The seal nosepiece must follow the seal seat surface 
under different operating conditions without surface contact or excessive 
increase in film thickness, which would yield high leakage. In addition, 
the self-acting lift pads give the seal a high radial stiffness enabling 
the seal to accommodate radial face deformations. An experimental investi- 
gation of self-acting lift pad performance, where face deformations occured, 
is reported in (6) . 

Since the seals must be pressure balanced, a proper balance of the 
opening forces (due to the pressure drop across the sealing dam and the 
force generated by the lift pads) and the closing forces (due to hydro- 
-static forces and spring forces) , must be achieved with a leakage gap that 
has tolerable mass leakage. The gap must be small enough so that the leak- 
age is minimal but it must be large enough so that power dissipation, due 
to shear in the film, and the inherent face deformations are tolerable. 

Thus, the design of the sealing gap is vital to seal performance, and the 
pressure distribution in the gap and mass leakage through the gap must be 
analyzed. 

In this paper, only the sealing dam portion of the seal (Fig. 1) will 
be analyzed. The classical viscous, isothermal, subsonic, compressible 
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flow analysis for parallel sealing surfaces is well known (7). Refer- 
ence (8) analyzed the parallel film constant area hydrostatic case includ- 
ing the effects of fluid inertia, viscous friction and entrance losses. 
Subsonic and choked flow conditions can be predicted and analyzed for both 
laminar and turbulent flows. Result^ showed good agreement with experi- 
ment. This paper will extend the analysis of (8) to include seal face de- 
formation effects on the force balance and mass leakage. This analysis 
should be an aid in both the design of ordinary gas film seals and espe- 
cially gas film seals with self-acting lift pads. 

The objective of the paper is to present mathematical analyses of 
compressible fluid flow across shaft face seals with face deformations. 

An approximate integral method will be used to analyze two models. First 
a quasi-fully developed flow model will be formulated which can be used 
whenever fluid inertia is negligible (subsonic flow); Secondly, a vari- 
able area flow model, where fluid inertia is considered, is shown to be 
valid for both subsonic and seal exit choked flow. Entrance losses are 
accounted for and turbulent flow can be analyzed by utilizing an appropri- 
ate friction factor-Reynolds number relation. 

Causes of Seal Face Deformations 

Distortions of the primary sealing faces are inherently present in 
gas film face seals. Distortions present include radial and axial dis- 
placements due to the centrifugal force, and are especially important 
under high rotational speeds as anticipated for advanced aircraft opera- 
tion. A typical centrifugal deformation is shown in Fig.. 2. Another 
common face deformation is thermal coning caused by an axial thermal grad- 
ient along the shaft. The hotter end of the shaft causes a differential 
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shaft radial displacement which results in the face coning illustrated in 
Fig. 3. Other distortions could be caused by: 1) pressure - due to high 
pressure drops and improper seal balance diameter; 2) mechanical; 3) asym- 
metry of rotating seal seat; and 4) tolerance buildup due to fabrication- 
and assembly. Generally, for internally pressurized seals, the distortions 
will cause divergent seal faces. 

Basic Model 

The sealing dam model (Fig. 4) consists of two coaxial circular rings 
separated by a very narrow gap. The sealing surfaces are radially deformed. 
A pressure differential exists between the rings' inner and outer radii. 

The fluid velocities are small in both the inner-diameter cavity and outer- 
diameter cavity which bound the sealing dam. 

It is essential to have very flat and parallel surfaces for satisfac- 
tory gas film seal operation (e.g., specified flatness of the sealing sur- 
face within two light bands of sodium). This is necessary to minimize oper- 
ating distortions. Heat transfer analyses and subsequent stress analyses 
indicate that relative face deformations of less than two milliradians can 
be expected (see (5)). Hence, the analysis is representative of small face 
deformations or tilts. Although the face deformation is never in reality 
strictly linear, a linear deformation will be assumed since a closed form 
explicit solution can be readily obtained in some special cases. In prac- 
tice, no two designs would ever have the same deformation anyway. However, 
the effect of relative surface deformation can be represented by an "effec- 
tive" or "apparent" linear tilt of the surface. 

As shown in Fig. 4, the effective tilt or deformation of the surface 
can be represented by a relative tilt angle, a, or by specifying the entrance 
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and exit film thickness. As shown in Fig. 4, for small linear tilts of 
the sealing faces, the film thickness at any distance along the leakage 
flow path X can be found from 


where 


h = h^ + ax 



tl] 

[ 2 ] 


Analysis 

Two models will be examined and formulated which will predict seal be- 
havior where an exact differential analysis model (9) is impossible to solve 
or impractical for design analysis purposes. Such cases are turbulent flow, 
where exact physical knowledge is unknown, and near-choked and choked flow 
where nonlinear behavior characterizes the flow. Approximate integrated 
average methods will be used. Although the integral models only satisfy 
mean conditions in the flow field, they may give good results on gross quan- 
tities such as seal leakage and pressure distribution. The constant area 
analysis (8) used an approximate integrated average method and showed good 
agreement with experiment. 

First a. quasi-fully developed compressible flow model will be presented 
for both parallel and small linear deformed surfaces. This model yields 
tractable solutions in relatively simple forms. Then a variable area analy- 
sis will be formulated for radial area expansion and small surface deforma- 
tions. This analysis will also include both the quasi-fully developed flow 
and constant area flow as special cases but the solutions will require a 
numerical Runge-Kutta solution. (The constant area analysis (8) used a 
linear iteration solution scheme.) 
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Quasi-fully Developed Flow Model 

This is the classical fluid flow case. This model is widely used to 
describe pipe and duct flows. This model is valid whenever the viscous 
forces dominate. That is, when entrance effects are negligible and the 
flow is subsonic (not near-choked or choked) . This flow can be called 
quasi-fully developed and is presented here from a seal point of view. 

Consider the control volume shown in Fig. 5 for situations when the 
fluid inertia is negligible. The momentum conservation is a balance be- 
tween the pressure and viscous friction force which is 

AdP = - t dA [3] 

w w 

Now introduce the following parameters 

hydraulic diameter, D = — — 

w 

dx 

_ x 

mean Fanning friction factor, f = 

pu 

2 


into Eq. [3] this results in 


D ~ = - 2pu 2 f 
dx 

Substituting the perfect gas law and mass flow definition 

M = pu A 

yields the following useful form 


[4] 


[5] 


P dP = 


-2f 1R T M dx 
DA 2 


[ 6 ] 
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Constant Area Flows 


Before analyzing sealing surfaces with deformed faces, several constant 
area cases will be solved. Assuming isothermal flow, constant area, mean 
friction factor, and hydraulic diameter, Eq. [6] can be readily integrated. 
The result is 


DA 


M = 


k- 


4f (RT(R 2 - R 2 ) 


[7] 


For radial flow between co-axial parallel disks and parallel plates, the 
hydraulic diameter, D, is given by 

D = 2h [8] 

Generally, the mean friction factor is related to Reynolds number by a rela- 
tion of the following form 


f = 


n 


[9] 


Re 

It is useful to express the Reynolds number in the following form 


„ 2M 
Re = — 
Wy 


[ 10 ] 


Now, both laminar and turbulent flow cases will be considered. 

1) Laminar Flow 

For laminar flow, the friction factor is derived from the classical, 
viscous compressible flow solution (9) and the derivation is also shown in 
(9) . The resulting mean friction factor - Reynolds number relation is 


f -i*. 

Re 


[ 11 ] 


Using relation Eq. [10] yields the following form for Eq. [6] 

„ , -12y (RT M , 

P dp = r - _ — dx 

Wh 


[ 12 ] 
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This equation is identical to a form which can be derived from the differ- 
ential analysis (9) . 

2) Turbulent Flow 

The Blasius relation of friction factor - Reynolds number appears to 
satisfactorily describe a large class of fully developed flows. Thus, in 
Eq. [9], k = 0.079 and n = 0.25. Substitution in Eq. [7] yields 


M = 


3.169 Wh 12/7 (p 2 - P 2 ) 

tK 4 / V/V / 7 (R 2 - V 4/? 


[13] 


which gives the functional relation of the variables in quasi-fully devel- 
oped turbulent flow. 

The pressure distribution can be found by integrating Eq. [6] from the 
entrance to any distance downstream, x. The result is 


P 

x 



0.1329 P 1/4 )RT M 7/4 x 

W 7/ V 


[14] 


This equation can be further simplified by substitution of the mass 
flow Eq. [13]. This yields 



[15] 


By examining Eq. [15], we see that the pressure distribution equation 
is the same as the one for laminar flow (9) and independent of fluid prop- 
erties. This suggests that it may be desirable from a leakage point of 
view to operate (if possible) in the turbulent flow regime. (For the same 
gap the pressure distribution is the same but the leakage is less.) 
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Variable Area Flows 


Eq. [6] can be integrated for both radial flow and constant width flow 
with small tilts of the sealing surfaces. Using the same restrictions as 
previously stated the following results are obtained for mass flow: 

1) Radial Flow, W = 3r 
(a) Laminar flow 


M = 


6h3 ( P 2 - P 1 

24pfRT In R 1 /R 2 


[16] 


Note, that if 3 = 2 tt, Eq. [16] becomes the same as the classical 
viscous, compressible, radial flow leakage equation found in (9). 
(b) Turbulent flow 


M = 


3.995 3 h 12/7 (p 2 


- P 


4/7 


pl/7 IR 4/7 x 4/7/ 1 


4/7 


3/4 t,3/ 4 , 


2) Flows With Small Linear Tilts and Constant Width, 
h = h^ + ax 
(a) Laminar flow 


[17] 



Wh ch»r( 1> 
24y IRT(R 


2 

1 

2 



[18] 


This equation is identical to the compressible viscous flow solution 

obtained from the differential analysis (9). The characteristic film 

thickness, h , , is defined as 

char 


h 


char 



[19] 
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Note when there are no tilts present h^ = h^ (hence a = 0) , then 

h , = h and the classical parallel surface mass leakage equation is 

char 

evolved. However, if the effect of an effective or linear apparent tilt 
is desired, the effect on mass leakage can be easily calculated from the 
parallel leakage equation. A simple computation enables leakage with de- 
formation to be readily calculated or read from a parallel leakage plot. 

For this case the radial pressure distribution across the sealing 
dam is 



Where the linear tilt factor, T 1 , is defined as 

h^(2h. + ax) 

T=— 7 [21] 

2h (h 1 + ax) 

m I 


Note for parallel surfaces a = 0, hence h = h = h 0 and the tilt 

m 1 2. 

factor, T = 1 which yields the parallel surface pressure distribution, 
(b) Turbulent Flow 


3.169 Wh 


i2/7/-2 _ ,2\ 


4/7 


M = 


char \ 1 


2 / 


V'V'V' 7 ^ - v 477 


[ 22 ] 


This Eq. [22] is identical to Eq. [13] for constant area and parallel 
surfaces except the characteristic film thickness again describes the small 
linear tilts. 

It is interesting to note that pressure profiles found with small face 
deformations present are still independent of fluid properties. However, 
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as will be seen in the next section, this will not be the case when convec- 
tive inertia effects are important. 

Variable Area Flow With Inertia 

The objective of this section is to present a mathematical analysis that 
includes fluid inertia, viscous friction, entrance losses, subsonic and 
choked flow conditions with area changes due to both radius change and/or 
small tilts. For subsonic viscous flows the analysis in the previous section 
can be used; however, with fluid inertia effects the following analysis must 
be used. The analysis will parallel that for the constant area flow in (8); 
however, it will be seen that resulting equation to be solved is more com- 
plex. The resulting friction parameter, Mach number, and area change equa- 
tion must be solved using a numerical solution scheme. The numerical scheme 
used here is the Runge-Kutta technique. 

The analysis can be separated into two parts, which can then be con- 
sidered separately. One part is an analysis of the entrance flow, while the 
other part is an analysis of the seal leakage path itself. 

Entrance Flow 

The entrance flow is treated identically as it was for constant area 
flow (8). That is, the entrance conditions are considered either isentropic 
or modified to account for entrance losses by an empirically determined 
velocity loss coefficient, C^. For example, the entrance pressure can be 
found from 
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Further details can be found in (8). For subsonic flows where fluid inertia 
is not important, such as the quasi-fully developed flow, the entrance con- 
ditions were shown to be negligible (8) (here Pq = P^) . However, entrance 
effects are important for choked flow conditions. 

Seal Leakage Passage Flow 

It is assumed that the flow in the seal leakage flow region behaves as 
a variable area adiabatic flow with friction. A quasi-one-dimensional ap- 
proximation is made wherein it is assumed that the flow properties can be 
described in terms of their cross-sectional averages. 

The following assumptions have been made in the analysis : 1) The ef- 
fects of rotation are neglected. 2) The flow is adiabatic. 3) No shaft 
work is done on or by the system. 4) No potential energy gradient is pres- 
ent such as caused by elevation differences, etc. 5) The fluid behaves as 
a perfect gas. 

The control volume is shown in Fig. 5. The governing equations with 
area changes reduce to the following differential forms: 

Conservation of Mass 


dp_ 1 du 2 

p 2 2 + 

H u 


dA 

A 


0 


[24] 


Conservation of Energy 

dl + <Y - 1) m 2 duf 

T 2 2 

u 

Equation of State 


0 


[25] 


dP = d£ dT 

P P T 

Conservation of Momentum (for a small area change) 


[26] 


-AdP - x dA = Mdu 
w w 


[ 27 ] 
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Introducting the hydraulic diameter and Fanning friction factor into 
Eq. [27] and combining Eqs. [24], [25], [26], and [27] results in a single 
equation to be solved which is 


dM 2 

-2 

r i 2 1 

1 + j(Y - DM J 

dA , 4fdx 

2 r i 2 

YM Z 1 1 + i(y - 1)M _ 


M 2 = 


1 - M 2 

" A + D 

1 - M 2 




[28] 


This is the identical result obtained from the Table of Influence Coeffi- 
cients for generalized one-dimensional flow in (11) and (12) . 

The dependent variables can be found by integrating directly from 
M to M* = 1, (A to A*, etc.) since the variables are separated. De- 
tails can be found in (9) and (12) . Performing the integrations results 
in the following equations: (Here the "starred" quantities denote the 

critical flow conditions.) 


£ i / 1 + 

“ " M \j i(Y ♦ 1) 


[29] 


T* 

~ 


1 + j(Y - DM 2 

|(Y + 1) 


[30] 



P* 

P 


am [ZEElEZ. 

A* V 1 + y(Y - 1)M 2 


A*M \j 1 + j(y - 1)M 2 


[31] 


[32] 


Solution Scheme 

The solution scheme is similar to the constant area case in (13) and 
will only be briefly described here. Basically the only known boundary 
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condition is that the flow is critical at some critical film thickness h . 

JL 

Hence, h is first determined. Then for h < h , the flow is assumed to be 

subsonic and for h > h*, the flow is assumed to be choked. The values of 

h and h considered are the characteristic film thickness, h , . The 

char 

entrance flow must match the seal passage flow. Hence, an iteration scheme 
must be used. The details of this procedure can be found in (13). 

Equation [28] must be solved numerically for the variable area cases 
of interest. The four point Runge-Kutta scheme is employed. Equation [28] 
is a first order differential equation in x. Therefore, one boundary con- 
dition on is required. The only known boundary condition is at x = x*, 

2 

M = 1. Hence, a fictitious length must be used for subsonic flow. It would 

be ideal to start at x and march backward to x = 0; unfortunately, the 
2 * 

derivative dM /dx must be known at x . Since the derivative is infinite 

£ 

at x , this location cannot be used as the starting point. The entrance at 
x = 0 is used as the starting point. An is guessed which will give 

M* - 1 at x*. 

The interval of integration is divided into subintervals. Since the 
Mach number is known to greatly vary in the last 10 percent of the seal 
prior to the choking point, the last subinterval is divided into smaller 
subintervals. Further details of this method can be found in (9). 

Results and Discussion 

Since experimental results are not known to exist for deformed seal 
surface flow under near-choked and choked flow conditions, some results 
will now be presented which were obtained from the variable area analysis. 
Both small linear tilts and pure radial area expansion flow will be consid- 
ered. (The analysis can be used for monotonically increasing or decreasing 
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film thickness distributions which would be a specified input.) All results 
presented here used isentropic entrance conditions and were in the laminar 
flow regime. The mean friction factor-Reynolds number relation given by 
Eq. [11] was used in the cases studied. The solutions were obtained using 
the Lewis Research Center IBM 7094/7040 direct couple computing. system. 

The approximate integral analysis was first used to solve a class of 
small surface deformation problems where the viscous differential model (9) 
was valid. The approximate integral and differential analysis showed good 
agreement. 

Figure 6 compares the pressure profiles obtained from using this anal- 
ysis with the viscous flow solution using Eq. [20] for a convergent seal 
face deformation. The conditions used in this analysis are representative 

of aircraft engine idle and are the same conditions used in the design ex- 

2 

ample in (10). The conditions, were : P^ = 65 psia (45 N/cm abs) , = 15 

psia (10.3 N/cm 2 abs), T Q = 100 F (311 K) , R =» 2.315 in (5.880 cm), and 
R 2 = 3.315 in (8.410 cm). The mean film thickness is 0.3 mils (7.6 ym) . 
Notice that there is excellent agreement between the two solutions along 
the first 40 percent of the seal leakage path (represented by a radial dis- 
tance from 0 to 0.020 in). Deviation from the small tilt analysis increases 
as the exit (0.050 in radial distance) is approached (the exit pressure, 
however, is held fixed as a boundary condition). This result is expected 
physically. In the first part of the seal leakage passage the flow is 
still primarily viscous; however, the flow tends to accelerate as it flows 
towards the exit. Hence, the fluid inertia effects become a concern and 
the two solutions are not expected to agree in this region. 
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Figure 7 again shows pressure distribution results obtained from the 
variable area approximate analysis for a divergent tilt of 2 milliradians . 
Distributions for mean film thicknesses of 0.1 mil (2.5 ym) , 0.2 mil (5.1 ym) , 
0.3 mil (7.6 ym) , and 0.5 mil (12.7 ym) are presented. The other conditions 
were Pq = 215 psia (148 N/cm^ abs), = 15 psia (10.3 N/cm^ abs), = 800 F 
(700 K) , R a 3.265 inch (8.300 cm), R^ = 3.315 inch (8.410 cm). These con- 
ditions are representative of advanced aircraft cruise conditions (10). These 
conditions represent subcritical (subsonic), critical (Y ^ = P^ and = 1) 
and supercritical (choked) flow conditions. Also shown is the parallel film 
pressure profile for 0.1 mil (2.5 ym) film thickness. This is the classical 
parabolic profile for viscous compressible flow. In addition. Fig. 7 shows 
a supercritical flow pressure profile for parallel sealing surfaces and a 
film thickness of 0.5 mil (12.7 ym) which. was obtained using the constant 
area analysis of (8) . The variable area approximate analysis shows excell- 
ent agreement with this parallel film profile with a 0.5 mil (12.7 ym) film 
thickness . 

Figure 8 shows the pressure distribution results obtained from the vari- 
able area approximate analysis under identical conditions as those in Fig. 7 
except a convergent tilt of 2 milliradians is considered. The parallel film 
pressure profiles are again presented as a reference. For critical flow the 
convergent film pressure profile indicates a substantially higher opening 
force than the parallel film profile. Notice that the variable area analysis 
again agrees with the constant area analysis case for the supercritical fl ox- 
case and a 0.50 mil (12.7 ym) parallel film. 

Other values of interest in gas film seal design are opening force, 
center of pressure and leakage flow rates. Table 1 shows the opening force. 
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center of pressure and entrance and exit Mach numbers for these divergent, 
parallel and convergent sealing face conditions. These values are for an 
operating point representative of an advanced aircraft under cruise condi- 
tions. Table 2 shows a mass leakage rate comparison for all three sealing 
surface cases. Both the divergent and convergent seal surface cases have 
the same characteristic film thickness but the leakage rates are not iden- 
tical, as is the case when the viscous forces dominate (classical viscous 
compressible flow). The convergent film case always has less leakage than 
the divergent film. This is due to the fluid inertia effect which is more 
pronounced for the convergent film where the fluid accelerates more rapidly 
than the fluid decelerates in the divergent film case. Also shown in 
Table 2 is the parallel film case where the mean film thickness equals the 
characteristic film thickness. The 0.500 mil (12.7 pm) mean film thickness 
case, only, is close to the tilt cases’ characteristic film thickness of 
0.497 mils (12.6 pm). Here the leakage for the parallel film is about half 
the difference of the convergent and divergent seal surface cases. 

Figure 9 shows the pressure distribution for pure radial viscous flow 
found from this variable area approximate analysis and the solution using 
Eq. [16]. The conditions were representative of aircraft idle operation: 

Pq = 65 psia (45 N/cm^ abs), = 15 psia (10.3 N/cm^ abs), T^ = 100 F 

(311 K) , R^ = 2.315 inch (5.880 cm), and = 3.315 inch (8.410 cm). The 
parallel surface case of 0.50 mil (12.7 pm) film thickness was solved. 

The variable area approximate analysis slightly underestimates the pressure 
along the seal passage length. This slight discrepancy is probably due to 
the choice of mean friction factor. The mean friction factor used here was 
the same as used for the constant area cases but referenced at the mean 
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radius. A problem arises in the selection of a mean friction factor for 
variable area flows. In pure radial flow the mean friction factor varies 
with radius (since the Reynolds number varies with radius) . The proper 
choice of this friction factor has to be examined in further detail exper- 
imentally. The friction factor - Reynolds number relation selection will 
also be a problem for large surface tilts and should also be investigated 
experimentally. 

As stated in (8), in most face seals the area expansion is negligible 
due to the radius ratios being close to one. Hence, the flow will usually 
choke at the exit. However, choking can occur at entrance. If there is a 
separation bubble present at the entrance (due, e.g., to the flow turning 
into the sealing faces) , there could be a large entrance area decrease with 
choking occurring at the vena contracta. In this case the analysis pre- 
sented here cannot be used in its present form. For smaller radius ratios 
such as those that characterize externally pressurized gas bearings, chok- 
ing will also occur at the entrance. 

SUMMARY OF RESULTS 

Approximate integral analyses have been presented for compressible 
fluid flow across shaft face seals with small face deformations. The fluid 
properties are averaged across the fluid film. This quasi-one-dimensional 
integral analysis includes fluid inertia and entrance losses in addition to 
viscous friction which is accounted for by a mean friction factor. Subsonic 
and choked flow conditions can be predicted and analyzed. The model is valid 
for both laminar and turbulent flows. The following pertinent results were 
found : 
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Quasi-fully Developed Flow 

a) Results for parallel films and small linear tilts agree with the 
exact classical viscous compressible flow solutions including pure radial 
flow. 

b) Leakage formulas were developed for quasi-fully developed turbulent 

flow. 

c) Small linear tilts of the sealing surfaces can be described using 
the parallel film leakage equation provided the characteristic film thick- 
ness is used. This is true for both laminar and turbulent flows. 

d) Pressure profiles are independent of fluid properties for both 
parallel surfaces and surfaces with small linear tilts. This is true for 
both laminar and turbulent flow. 

Variable Area Flow with Inertia 

a) Results agree with pure radial flow results and constant area flows 
with small seal face deformations. 

b) Fluid inertia affects both the pressure distribution and mass leak- 
age flow rates for divergent, convergent, and parallel sealing surfaces. 
Hence, for severe operating conditions, such as the advanced aircraft cruise 
conditions considered in this paper, inertia effects must be accounted for 
to get predictable gas film seal performance. 
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Table 1 - Comparison of various physical quantity values for divergent (+0.002 radian tilt), parallel, and convergent (-0.0Q2 radian 
tilt) sealing faces representative of subcritical, critical, and supercritical flow; P 0 ° 215 p.s.i.a. (148 N/cm 2 abs) , = 15 

p.s.i.a. (10.3 N/cm^ abs), Tq q 800° F (700 K) , Ri = 3.265 in. (8.300 cm), R2 0 3.315 in. (8.410 cm). These conditions are repre- 




Divergent tilt (+0.002 radians) 0.0068(0.0031) 0.0837(0.0380) 0.2800(0.1270) 0.9720(0.4410) 
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Fig. 5 - Model and notation of sealing faces including control volume for quasi 
one-dimensional flow with area change. 
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PRESSURE, P, LBF/IN. 2 -ABS 

Fig, 6 - Comparison of variable area approximate analysis with exact 
compressible viscous flow solution; negative one milliradian tilt, 

0.3 mil (7.6 pml mean film thickness, Pj ■ 65 psia (45 N/cm 2 absl, 

15 psia (10.3 N/ cm 2 abs), T 0 - 100° F (311 K), 3.265 in. (8.300 cm), 

and R 2 " 3.315 in. (8.410cm). 


o h » 0. 1 MIL (2. 5 pm) - SU8CRITICAL FLOW 

a h m - 0.2 MIL (5. 1 pm) - CRITICAL FLOW 

h m ■ 0.3 MIL (7.6 pm) - SUPERCRITICAL FLOW 

o h ■ 0. 5 MIL (12. 7 pm) - SUPERCRITICAL FLOW 

h m ■ 0. 1 MIL (2. 5 pm) - SUBCRITICAL FLOW I PARALLEL FILM CONSTANT 

h„ ■ 0. 5 MIL (12. 7 pm) - SUPERCRITICAL FLOWJ AREA ANALYSIS 

o h m ■ 0.5 MIL (12, 7 pm) - SUPERCRITICAL FLOW, PARALIEL FILM USING 

VARIABLE AREA ANALYSIS 
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Fig. 7 - Results using the variable area approximate analysis 
for pressure distributions; positive two milliradian tilt, condi- 
tions represent subcritical, critical, and supercritical flow, 
mean film thicknesses of 0. 1 mil (2.5 pm), 0.2 mil (5. lpm), 

0.3 mil (7.6 pm) and 0.5 mil (12.7 pm), Pg ■ 215 psia (148 N/cm 2 
abs), 15 psia (10. 3 N/cm 2 abs), Tg = 800° F 1700 K), Rj- 3.265 in. 
(8. 300 cm I, R 2 = 3.315 in. (8.4106 cm). Also shown are compa- 
rable parallel film cases. 
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o = 0. 2 MIL (5. 1 tim) - CRITICAL FLOW 
a h = 0.3 MIL <7.6 pm) - SUPERCRITICAL FLOW 

o h m = 0. 5 MIL (12. 7 pm) - SUPERCRITICAL ROW 

= 0. 2 MIL (5. 1 pm) ~ CRITICAL FLOW 1 PARALLEL FILM CONSTANT AREA AMAI YSIS 

h ffl » 0.5 MIL (12. 7 pm) - SUPERCRITICAL ROW ] HAKALLa HLM ' C0N5fANT AREA ANALYSIS 

o = 0.5 MIL (12.7 pm) - SUPERCRITICAL ROW, PARALLEL FILM, VARIABLE AREA ANALYSIS 



Fig. 8 - Results using variable area approximate analysis for 
pressure distributions; negative two milliradian tilt, conditions 
represent subcritical, critical, and supercritical flow, mean 
film thicknesses of 0.2 mil (5. lpm), 0.3 mil (7.6 pm), and 
0 . 5 mil (12 7pm), P 0 “ 215 psia (148 N/cm z -abs), P 3 » 15 psia 
(10.3N/cm 2 -abs), T 0 “ 800° F (700 K), Rj- 3. 265 In. (8.300 cm), 
R 2 “3. 315 in. (8.410 cm). Also shown are the parallel film 
cases. 
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Fig. 9 - Comparison of variable area approximate analysis 
with exact compressible viscous flow solution for pure 
radial flow; parallel film, 0.5 mil (12.7 pm) film thickness, 
Pq ■ 65 psia (45 N/cm 2 abs), P 3 = 15 psia (10. 3 N/cm 2 abs), 
Tq * 100° F (311 K), Rj * 2. 315 in. (5.880 cm), and R 2 « 

3. 315 in. (8,<tl0cm). 
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